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Abstract 

I present the theory and analysis behind the experiment by Fomalont 
and Kopeikin involving Jupiter and quasar J0842+1845 that purported 
to measure the speed of gravity. The computation of the vj/c correction 
to the gravitational time delay difference relevant to the experiment 
is derived, where vj is the speed of Jupiter as measured from Earth. 
Since the vj/c corrections are too small to have been measured in the 
Jupiter /quasar experiment, it is impossible that the speed of gravity 
was extracted from the data, and I explain what when wrong with the 
data analysis. Finally, mistakes are shown in papers by Fomalont and 
Kopeikin intended to rebut my work and the work of others. 
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1. Introduction 

Albert Einstein constructed his general theory of relativity so that classical gravity 
would be compatible with the principles of special relativity. As such, gravitational 
waves and the influences of gravity are suppose to propagate at the speed of light c. 
For example, if hypothetically the Sun were to explode into two pieces then the force 
of gravity on the Earth would change. However, it would not happen suddenly but 8 
and 1/3 minutes later, since this is the time it takes gravitational effects (and light) 
to travel from the Sun to the Earth. In other words, the Earth would continue in its 
almost circular orbit as governed by the gravity of a single massive central body for 
another 8 and 1/3 minutes, only after which would its motion be determined by the 
two exploding pieces. Since gravitational waves have not yet been detected, it has not 
been possible to test whether they travel at the speed c, nor has there been a system 
in which the speed of propagation of gravitational influences Cg has been measured. 

The lack of a measurement of Cg inspired S. Kopeikin to propose an experiment 
to test whether Cg is the speed of hght.[l] The ideas behind his proposal are outhned 
in the next four paragraphs. 

When electromagnetic waves pass by a massive object M, two effects occur: 
Firstly, the waves are very slighly bent, and, secondly, there is a tiny delay in the 
transmission time. Both are prominent effects of general relativity, and the latter is 
known as the Shapiro time delay. [2, 3, 4] 

The physical solution of Einstein's equations involve the position Xm of gravity- 
generating objects at retarded times: 

tret^t-\x-XM{tret)\/c. (l.l) 

Here, x is the location at which gravity is exerting its influence. The use of retarded 
times tret as given in Eq. (1.1) implies that the effect of gravity propagates at c. To 
allow for gravity to propagate at a different speed Cg, one would expect to replace 
Eq. (1.1) by 

tret — t ~ \x ~ X M{tret)\ / Cg ■ (l-^) 

On September 8, 2002, a conjunction of Jupiter and quasar J0842-I-1845 took 
place. Kopeikin argued that this event could be used to measure the speed of prop- 
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agation of gravity. [1] The thinking behind his proposal originates from the previous 
paragraph. Being very far away, the location of quasar J0842+1845 in the sky is 
virtually fixed. Jupiter, however, moves. Its position at the retarded time depends 
sensitively on its velocity and on Cg through Eq. (1.2). Therefore, one might expect 
that a precise measurement of the Shapiro time delay due to Jupiter on the signal 
from quasar J0842+1845 would permit a determination of Cg. 

The possible effect of Jupiter's velocity vj on the Shapiro time delay can be 
illustrated from an example: Suppose that Jupiter is moving toward the direction 
of the quasar waves. If Cg were infinite, which corresponds to the Newtonian limit 
of general relativity, then it would appear that the instantaneous position of Jupiter 
would be relevant. If Cg — c, then the position of Jupiter evaluated at the retarded 
time is somewhat farther away than the instantaneous position, and if Cg < c, then 
the position as determined by Eq. (1.2) is even further away. See Figure 1. 

Using an array of radio telescopes that stretched across the United States all the 
way to Germany, E. Fomalont and S. Kopcikin measured the tiny Shapiro time delays. 
By applying Very Long Baseline Interferometry (VLBI), they achieved a remarkable 
sensitivity at the picosecond level. The results for Cg, which were announced in Jan- 
uary 2003 meeting of the American Astronomical Society held in Seattle, Washington, 
immediately caught the attention of the media. The New York Times, for example, 
featured an article entitled "Einstein Was Right on Gravity's Speed." [5] At the meet- 
ing, Fomalont and Kopeikin announced that = c to within 20% and this result was 
subsequently published. [6] 

If this result is correct then it is a fundamental confirmation of Einstein's gen- 
eral theory of relativity. However, shortly after the American Astronomical Society 
meeting, a debate among astrophysicists arose about the Jupiter/quasar experiment. 
Several papers appeared arguing that Fomalont and Kopeikin had not accomplished 
their goal of measuring Cg. H. Asada had published an early work[7] stating that 
Kopeikin's idea actually measures the speed of hght instead of the speed of gravity, 
G. M. Will also argued that the measurements were not directly sensitive to Cg,[8\ and 
other papers also appeared criticizing the theory behind the experiment. [9, 10, 11] 
The difficulty that many of the above works were addressing is that there is no agreed 
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upon method for extending Einstein's theory to the case for which Cg ^ c. Kopeikin 
has repeatly tried to defend these criticisms of his work. [12, 13] 

In a Physical Review Letter [14] (hereafter referred to as my Physical Review 
Letter) , I bypassed the issue of how to extend general relativity to the case where the 
speed of gravity does not equal the speed of light and computed the Vj/c corrections to 
the Shapiro time delay in Einstein's theory using a relatively simple method. The Vj/c 
corrections did not agree with those derived by Kopeikin[l, 12] when his parameter 
Cg was set to c. The correct theoretical formula implied that the vj/c dependence 
was at least 100 times smaller than could have been measured by the array of radio 
telescopes used in the Fomalont /Kopeikin experiment. In other words, the speed of 
gravity could not have been extracted from the Jupiter/quasar measurement. My 
Physical Review Letter definitively settled the speed of gravity controversy: The 
parameter Cy has not been measured. 

Kopeikin had argued that there is an enhancement in the vj/c correction to the 
Shapiro time delay by a factor of 1/^^ where 9 is the angle between Jupiter and the 
quasar. [1, 12] Since 6 is small, it would seem to be that the vj/c correction is sizeable. 
However, no such 1/9 enhancement is actually present. 

My Physical Review Letter pinpointed the source of the discrepancy. The leading- 
order, velocity-independent part of the term that Fomalont and Kopeikin measured 
depends on the distance ^ of closest approach of the radios waves to Jupiter as 1/^. 
This distance is determined by the position of Jupiter and the radio waves as the 
latter pass by Jupiter. If Ri^j and 9obs respectively denote the Earth- Jupiter distance 
and the angle that an astronomer observes between Jupiter and the quasar, then 

C — REj9obs ■ (1-3) 

The angle 9obs is the one determined by the geometry of the positions of Jupiter, the 
quasar and the Earth at the time in which the radio waves pass by Jupiter. 

Fomalont and Kopeikin parametrized their data in terms of an angle 9i determined 
by the geometry of the positions of the above three objects at the time in which the 
radio waves arrived on Earth. These two angles differ: During the time in which the 
quasar signals travel from Jupiter to Earth, Jupiter moves a significant distance. The 
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relation between the two angles is 



Oobs ~ ^1 + 



n ■ vj 



(1.4) 



c 



where n is a unit vector pointing from Jupiter to the quasar's radio waves at the time of 
closest approach. When the factor 1/,^ appearing in the leading velocity-independent 
term is expressed in terms of 9i, it appears to become velocity dependent because 



In summary, parametrizing the Shapiro time delay using 6i makes the leading term 
seem to depend on the velocity of Jupiter. Furthermore, the fictitious velocity- 
dependent term appears enhanced. 

Fomalont and Kopeikin took their data and fit them to the leading order term but 
parametrized them in terms of 6i. They then extracted the vj/c dependence calling 
it the vj/cg correction. Given this procedure and that the data have error bars, it 
is not surprising that such a procedure produced the purported result that = c to 
within 20%. The conclusion that the speed of gravity is the speed of light to within 
experimental errors is not valid due to faulty analysis and a flaw in the theoretical 
understanding of the situation. 

Since the difference between Bobs and 9i is due to the change in the position of 
Jupiter as the quasar's radio waves propagate from Jupiter to Earth, it is clear that 
the parameter c in Eq. (1.5) is the speed of light and has nothing to do with the speed 
of gravity. Therefore there was no justification in using Cg in lieu of c when expressing 
the leading term in terms of 6i. 

The measurement of Fomalont and Kopeikin of the Shapiro time delay due to 
Jupiter is a remarkable experimental achievement. One should remember that the 
non-Newtonian effects of general relativity due to a planet had hitherto never been 
detected. However, the experiment has little theoretical significance or fundamental 
importance, and it indicates nothing about the speed of gravity. 

Most of our notation conforms to that of references [1], [12] and [14]. There are sev- 
eral small dimensionless parameters characterizing the Jupiter /quasar measurement: 
GnMj/{^c^) f» 6 X 10~^, where and Mj are respectively Newton's constant and 




(1.5) 
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the mass of Jupiter; vj/c 10~^; B/^ < 0.006 where B is the distance between 
any two VLBI stations on Earth; and Ogbs — ^/Rej ~ 0.001. These parameters, in 
the order given above, respectively represent the weakness of Jupiter's gravity on the 
radio waves, the non-relativistic nature of Jupiter's motion, the diminutive size of 
Earth compared to Jupiter, and the small observational angle between Jupiter and 
the quasar. 

In our analysis, we neglect the square of any of the above quantities. In particular, 
relativistic and high-order gravitational effects are ignored. 

2. The Leading Order Result 

Suppose that Jupiter is not moving. We refer to this as the static situation. Then, 
the Shapiro time delay for an electromagnetic wave travelling from the quasar past 
Jupiter to Earth is 



^^ = ^%^(l + l-(^^f^)) ' (2-1) 



where Rjq is the distance from Jupiter to the quasar. Eq. (2.1) is a textbook 
result. [15] 

In the quasar /Jupiter experiment, a scries of radio telescopes detected the quasar 
signals during the conjunction. The time difference A{ti,t2) between two such 
Shapiro delays At2 and A^i was measured: 

A(ti,t2) = At2- Ati. (2.2) 

The experimental situation is shown in Figure 2. Signals 1 and 2 propagate from 
the quasar past Jupiter and arrive at times ti and t2 on Earth at detectors located 

at positions Xi(ti) and X2(^2)- 

Using Eq. (2.1) in Eq. (2.2) yields 

A {t,M = At, - At, = In (^-^ J ^ ' ^^'^^ 

where A^ = — S ^-^d ru (respectively, r2j) is the distance between the first 
(respectively, second) detector and Jupiter. The last equality in Eq. (2.3) follows 
because A,^ is significantly smaller than either or ^2, and the differences in the 
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distances ru and r2j between Jupiter and detector can be neglected. We use ^ 
without a subscript to denote either of the detector-specific impact parameters or 
^2 when the distinction between the two is not important. Although gravitational 
effects are often long-ranged, the Shapiro time delay difference A (ti,t2) is generated 
in the vicinity of Jupiter as is evident from Eq. (2.3): A {ti, ^2) depends only on impact 
parameters. 

It is convenient to express A^ = Ci~ in terms of the displacement between the 
two detectors B — ^2(^2) — ^i(^i) because B is easily determined experimentally. See 
Figure 2. The radio signals from the quasar are bent slightly by an amount A(fi as 
they pass by Jupiter. However, it turns out that this effect can be neglected as we 
now show. 

The bending of a single wave is given by[15] 



A.^^^. (2.4) 



The angle that eventually arises between the two rays is 



SAip = Aip2 - Aip, = ^^^^^^^ . (2.5) 



Since the separation between the rays starts as A^ and increases as the distance times 

6A^, 



-n.B = AC + REjSA^ = AC{l + ^^^^^^]^AC. (2.6) 



The last equahty follows because 



AGnMjRej ^ AGnMjRej „ 

ec- - — ' 

where Rj is the radius of Jupiter. Since the angular deflection caused by Jupiter is 
so small, the separation between the two rays remains essentially constant. Indeed, 
including the second term in Eq. (2.6) in our analysis below only leads to corrections 
proportional to Newton's constant squared. 

Substituting Eq. (2.6) into (2.3), one obtains the following for the static situation 

'^GnMju ■ B , , 
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where we have used ^ = dobs^Ej- Eq. (2.7) is expressed in quantities measurable on 
Earth. 

3. The vj/c Corrections 

Now consider the case in which Jupiter is moving with a velocity vj with respect 
to the Earth. In principle, one should start with Einstein's equation 
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Ri^^ ^Zgi^'^R^ _::z^Ti'\ (3.1) 

in which it!^,^, the curvature tensor, is related to the stress-energy tensor T^^^ . 

For a static Jupiter, the stress-energy tensor only has a "00" component: T^^ — 
Pj{x), where pj is Jupiter's mass density. The metric g^i, in spherical-like coordinates 
about the center of Jupiter is given by the Schwarzschild solution 

c^dr' = (l - ^^^) c'dt' - (l - ^^^) ' dr' - - sin^ 9d^' , (3.2) 

a result that is only valid exterior to Jupiter. 

For a moving Jupiter, the stress-energy tensor has additional components T°*, 
which are proportional to Vj/c and generate vj/c corrections, and T*-', which are 
proportional to VjVj/c^ and may be neglected because Jupiter's speed is considerably 
less than the speed of light. The easiest way to obtain T^", R^^ and g^j,^, for the 
non-static case is to construct the Lorentz transformation that takes a non-moving 
Jupiter and sends it moving with velocity vj. This Lorentz tranformation is then 
applied to the tensors of the static case. Having obtained the metric for a non-static 
Jupiter, one would then need to compute the Shapiro time delay for this case. 

The above procedure for computing the vj/c corrections to the Shapiro time de- 
lay difference A(ti,t2) is, in principle, the one adopted by Kopeikin and is quite 
involved. [16] However, some simplifications occur because Jupiter's gravity is weak 
and vj/c'is small. One begins with the Newtonian approximation to Einstein's equa- 
tions and incorporates the Vj/c effects as perturbative corrections, which is the com- 
monly used post-Newtonian approximation. 

There is a simpler way to proceed, however, and that is to adopt a reference frame 
in which Jupiter is static. In such a frame, the Earth moves at a velocity ve given by 

ve = -vj , (3.3) 
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while Jupiter's velocity is zero. The formula in Eq. (2.7) for the static case is then 
valid. One only needs to incorporate the effects of having moving observation stations 
on Earth into the Shapiro time delay difference. This procedure is valid because (1) 
Einstein's theory is Lorentz invariant and (2) during the time in which the quasar 
rays propagate from Jupiter to the Earth, Jupiter moves almost in a straight line with 
constant speed. The propagation period is sufficiently short that the orbital motion 
of Jupiter around the Sun is not important. Since the same is true for the Earth, 
observers on both planets are essentially inertial during the time scales relevant to 
the experiment. 

— * 

Because the Earth is moving, the distance Bsf between points 1 and 2 as measured 

— * 

in this static- Jupiter frame is not equal to B as measured on Earth. In other words, 
if the first quasar signal arrives at Xi{ti) at time ti, then the Earth will move a short 
distance during the time in which it takes the second signal to arrive at X2{t2)- Place 
two observers in the static-frame (meaning that they are not moving with respect to 
Jupiter) so that one is located at the point 1 at time ti and the another is at the 
point 2 at time ^2- Then use these observers to make the time measurements. Since 
the situation is completely static, the formula for the static case may be used. 
The difference between the times at which the two measurements are made is 

t2-h^ \x2{t2) - Xo\/c - \xi{h) - Xo\/c + A {h, ^2) ■ (3.4) 

Here, Xq is the position of the quasar, and 1x2(^2) — ^o\/c— 1x1(^1) —Xq\/c is the time 
difference that occurs when gravitational effects are absent. Eq. (3.4) is an alternative 
definition of A (ti, ^2)- 

The leading contribution to this time difference is 

t2-h^ + A{ti,t2), 3.5 

c 

where the first term is, in general, larger than the second and arises from the first two 
terms in Eq. (3.4). Here, K, which is perpendicular to n, is a unit vector pointing 
in the direction of the quasar as seen from Earth. Since during the time t2 — ti, the 
Earth moves a distance ve {h — ti), the displacement between detectors in the static 
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frame Bgf is not the same as that in the Jupiter-moving frame B but the two are 
related by 

K B 

Bsf = B + VE {t2 -ti)^B —VE + A {tiM) VE . (3.6) 

The Earth's motion gives rise to three vj/c effects. The first occurs because Bsf 
needs to be used in Eq. (2.7). When this substitution is performed, a correction of 
AGmMjUs} ■ veK ■ B/{^c^) is generated. The effect of the last term in Eq. (3.6) may 
dropped since it is proportional to G%. 

The second vj/c correction occurs if the Earth has any motion toward (or away 
from) Jupiter. For example, if the Earth were moving toward Jupiter, then station 
2 would be moving toward the quasar signal 2 during the time after station 1 had 
detected signal 1 but before station 2 had received signal 2. Station 2 would then 
record a smaller time delay than if the Earth had not been moving. In other words, 
the time delay is reduced (or increased) by an amount 6 A (ti, that is equal to the 

— * 

time it takes light to travel the distance determined by the difference between Bgf and 
B. The corresponding correction due to the second term in Eq. (3.6) is independent 
of Newton's constant and is a contribution to the first part of Eq. (3.4) that involves 
detector distance differences. The third term in Eq. (3.6) leads to 

— * 

SA{t^,h)^-^^^A{t^,h) . (3.7) 
c 

To convert the result to the Jupiter-moving frame, one substitutes ve — —vj- 
Combining the above two effects with the leading term, one finds[14] 

AGnMj ^(^ , K-vj\ , K-BnsfVj \ 
A{hM)^ e^P^ I +^~)+ "c ) ■ ^ ^ ^ 

The above result is written in terms of quantities as measured by an observer on 
Earth with the exception of nsf, which gives rise to a third effect. Because the Earth 
is moving with respect to the static frame, the direction of the quasar as observed in 

— * — * 

the two frames differ: Kgf ^ K + {n • Vj/c) n. Since n is defined to be perpendicular 
to K, it too differs in the two frames: Hsf ~ n—{n ■ vj/c) K. When n^/ is substituted 
into Eq. (3.8), the last term is cancelled: 

AGnMj ^ ^(^ , K-vj \ 
A (ii, t2) = — n --Bin — I . (3.9) 
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As indicated by a lack of subscripts sf, all quantities are now those measured by an 
observer on Earth. Equation (3.9) has no l/^o^s term so that there is no enhancement 
of the velocity-dependent effects. 

The leading term for A (^1,^2) in Eq. (3.9) is of order 10~^° seconds, well within 
the measuring capability of the Jupiter /quasar experiment. However, the Vj/c term 
is not bigger than 10~^^ seconds, which is more than 100 times smaller than what 
was measurable by the VLBI stations. Therefore reference [6] was insensitive to the 
vj/c term. Furthermore, this vj/c correction is masked by larger corrections that are 
suppressed by factors such as B/^ and 9obs compared to the leading term. The above 
analysis leaves no doubt that Fomalont and Kopeikin did not measure the speed of 
gravity. 

Here is an example of a correction down by the order of Bobs that can be up to five 
times bigger than the vj / c term. If B has a component in the direction of the quasar, 
then the differences in the distances ru and r2j in Eq. (2.3) generate the following 
correction: 

(3.10) 

c3 Rej ^ ' 

It is smaller than the leading term by a factor of 0.5^/ Rej ~ 0.59obs- 

4. The Kopeikin Formula for A (ti, ^2) 

Using the post-Newtonian approximation, Kopeikin obtained the following result 
for A(ti,t2)[12] 

A(t„t2)=(^l + ^j^^ln 

where fij{si) = xi — Xj{si) and r2j('S2) = X2 — Xj{s2) are respectively the distance 
vectors between the observation points 1 and 2 and Jupiter evaluated at the retarded 
times 

s^ = ti- \xi - xj{si)\/c, 

S2^t2-\X2-Xj{s2)\/C. (4.2) 

Note that Kopeikin generally works in a frame in which Jupiter is moving and the 
Earth is not so that Jupiter's position xj{t) varies with time t, while that of a detector 
station does not. 



ru (.si) + A' • ru jsi) 
r2j {S2) + K ■f2j {S2) 



(4.1) 
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The term in Eq. (4.1) involving K ■ vj/c is discarded by Fomalont and Kopeikin 
due to its smallness. It is the main vj/c correction in Eq. (3.9). 

At time s, let ^i(s) (respectively, ^2(5)) be the angle between rij(s) (respectively, 
rij(s)) and K (the unit vector pointing toward the quasar). Then 

ru {s) + K- ru (s) = riX(s)/2 + 0{9t{s)) , 

r2j {s) + K- f2j (s) = r2jel{s)/2 + 0{9l{s)) . (4.3) 

The times Si and S2 at which one is to evaluate rij(si) and r2j(s2) in Eq. (4.3) are 
considerably earher than the observation times ti and t2. Indeed, since \xi—Xj (si) |/c 
is about the time it takes a quasar ray to travel from Jupiter to Earth, si in Eq. (4.2) 
corresponds to when the quasar ray passes near Jupiter. The same is true for S2- 
Therefore, it is Jupiter's position at this moment that is revelant. This is physically 
reasonable since this is when the planet exerts its biggest influence on the rays. 

Combining the results of the previous paragraphs, one sees that Eq. (4.1) agrees 
with Eqs. (2.3) and (3.9) since 

r2jel{s)^e2/r2j. (4.4) 

5. The Faulty Analysis of the Jupiter/Quasar Experiment 

Since Eq. (4.1) is the leading order result in Eq. (2.7) up to unmeasurable cor- 
rections, and Eq. (2.7) has no vj/c dependence in it, the question arises as to how 
Fomalont and Kopeikin exacted Cg from their data. Although the details of the analy- 
sis have not been revealed, I have surmised what transpired through their publications 
and through correspondence with Sergei Kopeikin. 

Fomalont and Kopeikin express the result for the Shapiro delay time difference in 
Eq. (4.1) of 

rij isi) + K -fij (si) 



S (si,S2) = 5 — In 



r2j (S2) + K ■f2j (S2) 



(5.1) 
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which is a function of retarded times Si and S2, as a fuction of the observation time 
ti. The leading order contribution to -S" (si, S2) is 



In 



6(S2) 



(5.2) 



If Jupiter is moving toward (or away from) the quasar then the distances ^1 (t) and 
^2 (t) between Jupiter and quasar-ray-trajectories decrease (or increase) with time t 
and are smaller (or larger) if evaluated at the time of observation: 

Rej 



6 (si) ~ 6 (^i) +n-vj 
6 (S2) ~ 6 (^i) +n-vj 



c 

Rej 



(5.3) 



See Figure 2. 

As is physically clear, the difference A,^ = ^i{t) — ^2{t) docs not change with time 
as Eq. (5.3) shows. However, the leading order result in Eq. (5.2) now does since the 
substitution in Eq. (5.3) leads to 



In 



In 



.6 (S2) 

The second term in Eq. (5.4) is 
fi ■ vjRej (I 1 



6 {h) 



+ 



n ■ VjRej 



(5.4) 



A^ (ti) n ■ VjRej _ n ■ B n ■ vj 
" c ~ ej{h) cRej 



(5.5) 



Summarizing, 



S{s^,S2)- S{h,ti) = A 



AGnMju ■ Bn-vj 



R 



(5.6) 



RejOI {ti) 

The right-hand side of the equation contains the artificially 1/9'^ enhancement claimed 
by Kopeikin.[l, 12] It arises because the position of Jupiter changes as the quasar 
signals travel from the Jupiter region to Earth. The distance from Jupiter of the 
(almost) linear path of the quasar ray when it is detected at a VLBl station is different 
from the distance when it passed by Jupiter. Fomalont and Kopeikin used tables to 
determine Jupiter's position at ti. Such information allows one to determine S {ti, ti) 
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analytically. So if one experimentally measures S {si, S2) and subtracts S {ti, ti), then 
one can fit the data to Ar in Eq. (5.6) to extract the artificially generated second term. 

Since S{si,S2) has no measureable vj/c dependence in it, there is implicit vj/c 
dependence in S{ti,ti). This is because 



AGnMju-B 



can be written as 



AGnMju-B 

S (ti.ti) ^ (5.8) 



to order (f j/c)^. 

The Introduction argued that perhaps one way to extend Einstein's theory to the 
Cg c case is to replace 1/c by 1/cg in the evalution of retarded times. Kopeikin does 
this in S (si, S2) when parametrizing it in terms of ti. This, however, is incorrect. The 
difference Ar between S{si,S2) and S {ti,ti) is due to Jupiter's change of position 
during the time it takes quasar rays to travel from the vicinity of Jupiter to Earth. 
This depends on the speed of rays, which is the speed of light, and not on the speed 
of gravity. 

Thus, when data with error bars are fit to Ar and the fitting function uses Cg 
instead of c, one is guaranteed to obtain the result Cg ~ c. Fomalont and Kopeikin's 
announcement that the speed of gravity is the speed of light to within 20% has no 
content. 

6. Can the Speed of Gravity Be Defined for the Jupiter/Quasar Exper- 
iment? 

This section addresses the theoretical issues raised in references [7, 8, 9, 10, 11]. 
Those references argued that Kopeikin's formula for the Shapiro time delay difference 
should involve the speed of light and not the speed of gravity. This debate is over 
how to define the speed of gravity in Einstein's general theory of relativity. 

In the static frame, Jupiter is not moving. The curvature of space-time created 
by the massive planet is governed by Eq. (3.2) and is static. The effects of gravity 
are not propagating and the speed of gravity concept is non-existent. The velocity 
dependent corrections obtained in Eq. (3.9) arise due to the motion of the VLBI 
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stations as they detect the radio waves. Since radio waves travel at the speed of hght, 
the velocity dependent corrections must be proportional to vj/c, and this is evident 
in the derivation of Eq. (3.9) in Section 3. It makes no sense to replace vj/chj vj/cg 
from the viewpoint of the static frame. 

On the other hand, in the frame in which Jupiter moves, the Introduction sug- 
gested that to extend Einstein's theory to the case Cg ^ c one should replace retarded 
times as determined by c by retarded times as determined by Cg. It would therefore 
seem as if the speed of gravity concept could be defined for this situation. 

However, results as measured in one frame must be consistent with those measured 
in another. One is forced to conclude that the speed of gravity concept is theoretically 
inconsistent for the Jupiter/quasar experiment. If one tries to define the general 
theory of relativity with Cg^ c for this problem, one violates Galilean invariance. 

Originally, C. Will has suggested that perhaps the speed of gravity might enter as a 
f j/c^ effect. [8] However, the above reasoning still applies. In relating the static frame 
results to the Jupiter-moving frame, one needs to use Lorcntz transformations if higher 
velocity effects are to be considered. Therefore, an extension of Einstein's theory to 
describe the Jupiter/quasar experiment for which the linear velocity corrections are 
of the form vj/c but the quadratic corrections are of order I'^/c^ for Cg ^ c would 
violate Lorentz invarince (but not Galilean invariance). Recently, C.Will has also 
come to the same conclusion that Cg does not appear in any higher power (vj/c)"' 
correction for a constantly moving Jupiter. [17] 

As mentioned in my Physical Review Letter, the above argument would fail if 
Jupiter (or another massive object) were accelerating toward (or away from) the 
quasar rays (or other electromagnetic waves). It is possible that the speed of grav- 
ity could be defined for this situation. The parameter Cg would then be associated 
with acceleration effects. It might be worth analyzing this case as a theoretical pos- 
sibility. Unfortunately, it is unlikely that such a system within or beyond the solar 
system exists with sufficiently large effects as to be measurable with current VLBI 
instruments. 
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7. The Response by Fomalont and Kopeikin 

In fairness to Fomalont and Kopeikin, it should be said that they have not ac- 
cepted the conclusions of my Physical Review Letter nor the criticisms of others. 
They have continued submitting papers[13, 18, 19] arguing that they did indeed mea- 
sure the speed of gravity. This section addresses those papers. Derivations of the 
formula in Eq. (4.1) mostly make up the content of these attempted rebuttals, while 
a few paragraphs are devoted to addressing the criticisms of other authors and of my 
Physical Review Letter. These paragraphs contain errors and false statements. 

For example, reference [13] says, "This part of the experiment was drastically mis- 
understood by Samuel who assumed that we measured position of quasar with respect 
to Jupiter by measuring the relative position of the quasar with respect to Jupiter in 
radio." Similar statements appear in reference [19] ("A fundamental flaw in Samuel's 
interpretation was his assumption that the direction to Jupiter was directly measured 
by VLBI network in the detection experiment so he confused the propagation of grav- 
ity and the propagation of radio waves." ) and reference [18] ( "Unfortunately, Samuel 
incorrectly assumed that the experiment directly compared the radio position of the 
quasar with that of Jupiter, and that the direction of Jupiter was determined by a 
photon reflected from its surface."). My Physical Review Letter never made such 
statements, nor does this review: The criticism of the manner in which the data anal- 
ysis was performed, which is presented in Section 5, focuses on the parametrization 
of and the expansion about the observation time ti by Kopeikin of the Shapiro time 
delay difference. 

The unpublished work [18] continues "The experiment monitored the position of 
the quasar as a function of the atomic time by the arrival of the quasar's photons 
at the telescope, while the Jupiter's position was determined separately via a precise 
JPL ephemeris, evaluated at the same atomic time as the arrival of a photon (via 
standard transformations from ephemeris time to atomic time). Hence the actual 
angle used for measuring A is 9i, not Samuel's 9obs- Thus, the vj/c correction 
was clearly separated from As and measured with a precision of 20%." Although 
these statements are intended to underpin my Physical Review Letter, anyone who 
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truly understands the situation will realize that they actually support the Letter and 
the criticism in Section 5 levelled at the data analysis of the Fomalont/Kopeikin 
measurement. 

Here is example of mistating the work of others as a means of defending the 
theory behind the quasar/ Jupiter experiment. Reference [19] says, "Our definition of 
the speed of gravity is more general than that used by Asada, Samuel, and Will [[7], 
[14], [8]] who hmited its meaning as the speed of propagation of gravitational waves." 
and "In their formulations of the experiment, these authors [Asada, Samuel, Will] 
assumed only far-field gravitational effects, where gravitational waves are dominant 
and differentiation between c and Cg occurs only at orders of {v/cf' beyond Shapiro 
delay and higher. This was one reason why the 'speed of fight' was interpreted as 
causing the observed aberration. However, the experiment was performed in the 
near-field of the quasar radio wave- Jupiter interaction where gravitational modes not 
associated with gravitational waves arc dominant." 

The work of the three abovc-citcd authors focused on the speed of propagation 
of the effects of gravity and not on the speed of gravity waves as falsely claimed in 
these quotes. My Physical Review Letter and this review did not discuss gravitational 
waves and have emphasized that the Shapiro time delay difference is due to relatively 
short- distance effects; I never "assumed only far- field gravitational effects" . In actual 
fact, it is Kopeikin's formalism that ends up distorting and mixing up the short- and 
long-distance gravitational effects of Jupiter as explained in Section 5. 

Reference [13] states, "The goal of the jovian defiection experiment was to distin- 
guish two angles 9i and 6obs- Confirmation that the apparent position of the quasar 
in the sky makes the angle 9i rather than 6obs with respect to Jupiter is a proof that 
gravity propagates with the speed Cg." This is not true: The difference between 9i 
and 6obs is due to the motion of Jupiter during the period in which the quasar waves 
travel to Earth. Quasar waves travel at the speed of light so that Cg cannot be ex- 
tracted from a difference of 9i and 9obs', See Eqs. (1.4) and (5.6) where it is incorrect 
to replace c by c^. 
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Reference [13] presents another example of faulty reasoning: "The equation 

does not contain terms being quadratic in 1/9. It may make an impression that 
the orbital motion of Jupiter does not provide any significant deviation from the 
Einstein's prediction of the light deflection because all velocity-dependent terms in 
the right side of Eq. (7.1) are smaller than the main term (proportional to {n ■ B)/6) 
by a factor of 10~^ and can not be observed with the present-day technology. This 
was the reason for Samuel's statement that terms of order v/c beyond the Shapiro 
time delay are not observable. This statement is erroneous because the Shapiro time 
delay must be calculated in terms of the present position of Jupiter at the time of 
observation ti." 

The last statement is incorrect. Kopeikin's own formalism (the right-hand side 
of Eq. (4.1)) leads to the conclusion that the Shapiro time delay difference A should 
be computed at the retarded time, which is near the time at which the quasar rays 
pass Jupiter. Furthermore, it is not necessary to evaluate A at ii: Since the Shapiro 
time delay difference is generated when the rays are in the vicinity of Jupiter, ti, 
which determines the retarded time si in Eq. (4.2), may be evaluated at any time 
after the rays have passed well beyond Jupiter. As is physically clear, tens of millions 
of kilometers beyond Jupiter, the Shapiro time delay difference will have been almost 
100% generated and then remain essentially unchanged. As is evident from Eqs. (2.3) 
and (3.10), there is only very weak dependence of A on the Earth- Jupiter distance 
Rej- 

In both references [13] and [19], Kopeikin claims that his formalism is Lorentz 
invariant, which is not the case. In Section 6, we argued that il Cg ^ c then not only 
is Lorentz invariance violated for the Jupiter/quasar experiment but also Galilean 
invariance. 

Kopeikin derived Eq. (7.1) in part to show that his formalism agreed with the result 
in my Physical Review Letter when Cg — c. Eq. (7.1) is Kopeikin's generalization to 
arbitrary values of the speed of gravity. In this equation, ve and vj are the velocities of 
the Earth and Jupiter as determined in a coordinate system in which the Sun is at rest. 
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However, Eq. (7.1) does not depend on vj—ve as required by Galilean invariance. The 
frame-dependence is manifest: Hypothetically place the Sun at another location, give 
it a velocity of vs relative to the original static Sun, have the Earth and Jupiter move 
with the same velocities that they originally had in the Jupiter/quasar experiment, 
then Eq. (7.1) becomes 

^GnMj ( Lj^ K -{vj- Vs) \ fi. B + K ■ B ( ^ ' (^^ ~ ^s) _ ft ■ {vj - vs) 
c^Orij yy Cg I V c Cg 

(7.2) 

because, with respect to the "new" Sun, Jupiter and the Earth now move with veloc- 
ities vj — Vs and ve — vs- The dependence on the "new" Sun's velocity vs does not 
drop out. How can the Shapiro time delay difference due to Jupiter and observed on 
Earth be so dependent on the Sun's motion? The answer is that Eq. (7.1) is wrong. 
In a frame in which both the Earth and Jupiter move, the correct result is Eq. (3.9) 
with vj replaced hj vj — ve- 

Reference [19] emphasizes incorrectly that it is the postion of Jupiter at the time of 
observation that is relevant: "the gravitational force, acting on photons as a space-like 
vector, is sensitive to the present position of Jupiter." For the sake of argument, let us 
assume that this is the case. When Kopeikin expands about the time of observation 
ti, he obtains, for the case Cg — c, avj/c correction of[l] 

Aih,h)^-^^(B-vj-K-vjK-B) . (7.3) 

It is easy to sec that this equation leads to physically unreasonable results. Suppose, 
for example, that Jupiter is moving toward the quasar rays. Hypothetically place 
a second planet twice as far away as Earth is to Jupiter and put VLBI stations 
on it along the same lines of observation as determined by quasar ray trajectories 
that pass through VLBI stations on Earth. Then one would expect to measure to 
a high precision the same Shapiro time delay difference. However, Eq. (7.3) predicts 
that more than twice the Earth-based correction will be measured. This follows 
because 1I{Q\t\j) — rij/^l{ti) where ^i(ti) is the distance between Jupiter and the 
worldline of the quasar ray 1 at the time of observation ti. In making the observations 
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"downstream" on the hypothetical planet, ru is doubled and ^i{ti) is smaller if 
Jupiter is moving toward the rays because the measurement is made later compared 
to the Earth-based one. How can Jupiter, which is so far away, have such a dramatic 
long-ranged influence on A? 

The answer is that vj/c corrections are not given by Eq. (7.3); The correct result 
is given in Eq. (3.9). Kopeikin should not have performed an expansion about the 
observation time that separated the leading Shapiro time difference result into two 
pieces. As explained in Section 5, his leading order result actually has implicit vj/c 
dependence in it (because it is written in terms of 9i instead of Bobs or which, 
when combined with his purported vj/c correction, gives the ^;j/c- independent result 
of Eq. (2.7) to order Vj/c^. This leading order term does not depend on how far 
"downstream" the measurement is made. The same is true of the vj/c correction in 
Eq. (3.9). 

8. Summary 

This work provides the leading vj/c corrections to the Shapiro time delay differ- 
ence and shows that they do not correspond to the ones used in the data analysis by 
Fomalont and Kopeikin even when the speed of gravity parameter Cg is set equal to c. 
The error made by Kopeikin is that he separated the leading f j/c-independent term 
into two pieces by expanding about the time of observation ti through the subtraction 
procedure involving Ar in Eq. (5.6). This introduced an artificially j/c-dependent 
term. It is enhanced by l/9i and is an artefact of the ti expansion. Although this 
vj I c-dependent term should depend the speed of light c because it is related to the 
change in position of Jupiter during the time in which radio waves travel to Earth, 
Kopeikin incorrectly replaces c by Cg. Hence, when data is parametrized in this way, 
one is guaranteed to obtain a result for Cg that is the speed of light to within ex- 
perimental errors. Finally, when Cg ^ c agree with others[7, 9, 10, 11] that the 
theoretical formalism[16, 21, 12, 20, 13, 18, 19] of Kopeikin and his co-workers is 
flawed when applied to Jupiter/quasar system and demonstrated this by showing 
that not only does it violate Lorentz invariance but also Galilean invariance. While 
the Jupiter/quasar measurement was an extraordinary experimental undertaking and 
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a fine achievement in precision, it had nothing to do with the determination of the 
speed of gravity. 
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Figure Captions 

Figure 1. The Expected Effect of the Speed of Gravity on the Gravitational Force 
Due to the Motion of Jupiter. 

Jupiter, shown as an open circle, is assumed to be moving toward the electromagnetic 
wave. In the top figure, Cg = oo and it is the instantaneous position of Jupiter at time 
t, which is given by the shaded circle, that is relevant. In the middle figure, Cg — c 
and the relevant distance is farther away: Since it takes time for the gravitational 
infiuence of Jupiter to propagate to the electromagnetic wave, it is Jupiter's position 
at an earlier time that is relevant. In the lower figure, Cg < c. 

Figure 2. The Definitions of Various Quantities Relevant in the Jupiter/Quasar Ex- 
periment. 

The diagram is not drawn to scale for reasons of clarity, and, in particular, angles are 
much larger than in the actual experiment. 
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